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Abstract 


Part  I: 

A  method  was  developed  to  find  the  far-field  radiation  pattern 
of  a  buried  horizontal  electric  dipole  (HED)  at  37.5  MHz.  The  imaginary 
part  of  the  index  of  refraction  was  shown  to  be  negligible  for  dry  soil 
at  this  frequency  so  standard  antenna  theory  and  ray-optic  theory  *e 
used.  The  effect  of  the  ground-air  interface  was  modeled  using  t 
transmission  coefficient  and  Snell's  law  for  a  dielectric  interf; 

Because  the  current  distribution  for  the  buried  HED  depends  on  anc  ua 
construction,  results  are  shown  for  the  far-field  pattern  in  the  air  for 
different  current  distributions  on  the  HED. 

The  literature  on  this  problem  was  reviewed;  most  used  the 
Sommerfeld  or  moment  methods  to  make  the  same  calculations.  The  results 
of  one  of  the  reports  using  the  Sommerfeld  method  could  be  compared  and 
were  found  to  be  similar.  An  extensive  bibliography  is  included. 

Part  II: 

The  analysis  was  then  applied  to  a  buried  antenna  array.  The 
current  distribution  was  known  and  was  used  to  calculate  the  far-field 
pattern.  It  was  concluded  that  the  far-field  pattern  is  highly  dependant 
on  the  current  distribution.  This  part  is  classified. 


I.  Introduction 


The  theory  of  antennas,  both  in  free  space  and  above  a  con¬ 
ducting  ground  plane  has  been  well  defined  and  tested.  When  an 
antenna  is  buried  in  a  medium  such  as  the  ground,  the  nature  of  the 
antenna  is  changed.  The  radiation  pattern,  current  distribution, 
impedance  and  other  parameters  can  no  longer  be  calculated  by  the  usual 
methods. 

Problem 

To  analyze  the  effectiveness  of  a  buried  antenna  array,  the  far- 
field  pattern  must  be  determined.  The  goal  of  this  thesis  is  to 
develop  a  simple  but  accurate  method  to  determine  this  pattern.  The 
antenna  array  to  be  studied  is  made  up  of  horizontal  electric  dipoles 
(HED's)  and  operates  at  37.5  MHz.  The  HED's  are  a  half-wavelength 
long  (4  meters).  The  exact  arrangement  of  the  array  is  classified 
and  is  discussed  in  Part  II. 

Definition  of  the  Far-Field 

The  far-field  condition  is  R>>D  where  R  is  the  observation 
distance  and  D  is  the  largest  dimension  of  the  antenna.  The  far- 
field  considered  in  this  study  does  not  include  the  fields  in  the 
ground.  Those  fields  are  attenuated  at  a  rate  greater  than  l/R.  The 
attenuation  is  dominated  by  the  lossy  effects  of  the  ground  over 
which  it  travels.  The  ground  wave  thus  is  insignificant  for  long  ranges 


y 

> 


^  it 


from  the  antenna.  Therefore  in  this  study  it  is  necessary  to  compute 
only  the  space  wave  for  the  far-field  pattern.  Its  fields  vary  as  a 
function  of  1/R  and  travel  away  from  the  surface. 


Development 

A  review  of  the  applicable  open  literature  is  given  in  Chapter  II. 
In  this  study,  the  radiation  equation  for  an  antenna  in  free  space  is 
applied  to  this  problem  using  ray-optics  and  Snell's  law  at  the  ground- 
air  interface.  This  is  presented  in  Chapter  III.  Chapter  IV  shows  the 
results  of  computations  of  the  pattern  for  different  parameters  using 
the  methods  from  Chapters  II  and  III.  Finally,  Chapter  V  gives  conclu¬ 
sions  and  recommendations  for  further  study. 
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II.  Literature  Review 


I 

There  are  many  studies  of  buried  antennas  in  the  open  literature. 
Many  of  them  analyze  the  buried  HED  at  lower  frequencies.  It  was  dif¬ 
ficult  to  apply  the  results  to  this  application  because  the  assumptions 
and  approximations  are  implicit  in  the  mathematical  derivations.  This 
is  primarily  due  to  the  complexity  of  modeling  the  conductive  half¬ 
space  of  the  ground.  Most  of  the  studies  reviewed  here  are  based  on 
Sommerfeld's  method. 

Sommer f eld  Method 

A.  Sommerfeld  did  the  initial  work  in  this  area  in  the  early 
1900's  (Ref  1)  and  his  method  is  followed  today.  A  quick  summary  is 
given  here  before  moving  on  to  the  specific  articles. 

The  Sommerfeld  method  obtains  an  exact  integral  representation 
of  the  fields  in  adjacent  conductive  half-space  and  lossless  half¬ 
space.  The  Helmholtz  equation 

(V2+k2)f  -  -jldln/k  (2-1) 

is  solved  by  a  Fourier-Bessel  transform  using  the  appropiate  boundary 
conditions.  The  fields  are  obtained  using 

E  -  V(V*f)  +  k2ff  (2-2) 

H  -  ' &  x  n  (2-3) 

n 

It  is  fairly  simple  to  show  that  equations  (2-1),  (2-2),  and  (2-3)  are 
a  solution  to  Maxwell's  equations.  The  inverse  Fourier-Bessel  transform 
can  then  be  found  by  an  asymptotic  approximation,  which  results  from  a 


suitable  integration  in  the  complex  plane. 

Several  reports  have  been  written  by  Biggs  and  Swarm  that  use 

the  Sommerfeld  method  (Ref  2-5).  They  calculate  both  the  near  and  far 

fields  from  a  buried  HED  operating  at  VLF  (3-30  kHz)  or  LF  (30-300  kHz). 

They  state  that  propagation  is  limited  to  low  frequencies  because  of 

losses,  but  it  is  not  expressed  as  a  formal  limitation  of  the  results. 

Biggs  (Ref  2)  analyzes  an  inclined  dipole  in  a  conducting  medium. 

He  assumes  that  the  depth  of  burial  is  much  less  than  the  observation 

distance.  He  also  assumes  that  n2»l  and  k  R/n2>l,  where  n  is  the 

g  a  g  g 

complex  index  of  refraction  of  the  ground,  k  is  the  propagation 
constant  of  the  air,  and  R  is  the  observation  distance. 

Biggs  and  Swarm  (Ref  3)  analyze  the  HED  in  a  conducting  medium. 
The  real  part  of  the  complex  index  of  refraction  is  neglected,  assuming 

-i-  »  e  (2-4) 

U)e0  r 

This  a  poor  approximation  at  VHF  (30-300  MHz)  as  will  be  shown  in 
Chapter  III.  Biggs  and  Swam  (Ref  4)  give  results  for  the  HED  in  a 
lossy  dielectric  consisting  of  ice  or  soil.  The  condition  in  this 
report  is  that 

n2 

R  sin0  »  r&-  (2-5) 

a  k 

a 

where  0^  is  the  angle  in  the  air  off  the  normal  to  the  ground.  This 
article  is  a  summary  of  an  earlier  report  (Ref  5)  that  goes  into  more 
detail  on  the  mathematical  analysis.  In  these  last  three  reports 
(Ref  3-5),  there  is  no  limitation  given  for  the  depth  of  burial. 

The  electric  fields  are  given  as  functions  of  polarization, 
distance,  and  position  for  all  4  reports  (Ref  2-5)  for  both  space  waves 


and  gro»'>'l  waves.  The  equations  are  similar  and  are  compared  in 
Chapter  IV. 

Entzminger  et_  al_  (Ref  6)  use  the  same  space  wave  equations  as 
given  in  Biggs  and  Swarm  (Ref  4.205)  to  calculate  fields  for  a  buried 
HED  operating  over  the  entire  HF  (3-30  MHz)  range.  They  also  study  the 
gain  and  impedance  of  an  insulated  wire  antenna.  Several  antennas  of 
this  type  were  constructed  and  buried  1  to  3  feet  deep  in  an  open  field. 
Both  the  space  wave  and  the  ground  wave  patterns  were  measured  for 
frequencies  of  2  to  10  MHz.  They  cite  excellent  agreement  between 
theoretical  and  measured  values. 

The  Sommerfeld  method  is  explained  in  detail  by  Banos  (Ref  7). 

He  calculates  both  the  near  and  far  fields  for  a  HED  at  any  depth  of 
burial.  He  assumes  LF  operation  in  a  conducting  medium  such  that  the 
real  part  of  the  complex  index  of  refraction  is  negligible.  (  A  special 
note  -  'n'  in  his  book  represents  the  reciprical  of  the  index  of 
refraction  for  the  ground.) 

Vaziri  (Ref  8)  uses  the  Sommerfeld  method  to  analyze  all  four 
Hertzian  dipoles,  horizontal  and  vertical  electric  dipoles  and 
horizontal  and  vertical  magnetic  dipoles,  as  buried  antennas.  Rather 
than  obtaining  the  inverse  Fourier-Bessel  transform  asymptotically, 
he  evaluates  the  inversion  integral  by  direct  numerical  integration  in 
the  complex  plane.  His  Fortran  program  is  given  but  its  limitations 
are  not  discussed.  Its  purpose  is  to  calculate  the  E  and  H-fields  at 
any  polarization  at  any  observation  point  in  a  certain  plane.  For  the 
HED,  this  plane  is  the  vertical  plane  parallel  to  and  intersecting  the 
dipole.  Above  ground  antennas  can  also  be  modeled  by  switching  the 
medium  constants  around. 


II-3 


Experimental  results  are  given  for  a  horizontal  traveling  wave 
antenna  buried  40  centimeters  deep  in  an  open  field  and  operated  at 
144  MHz.  The  theory  and  measurements  agreed  within  the  "acceptable 
accuracy".  Also,  Banos's  approximations  are  computed  in  the  program 
and  are  also  found  to  agree  with  the  measurements  in  most  cases.  This 
is  surprising  because  they  are  not  supposed  to  be  valid  for  such  high 
frequencies. 

King,  Sandler,  and  Shen  (Ref  9)  use  the  Sommerfeld  method  to 

find  the  fields  in  the  ground  for  a  buried  HED.  The  fields  in  the  air 

axe  not  examined.  Both  a  complex  permittivity  and  complex  conductivity 

are  used  to  calculate  the  propagation  constant  of  the  ground,  k  . 

8 

King  and  Shen  (Ref  10)  use  the  Sommerfeld  method  to  calculate 

the  fields  directly  above  a  buried  HED  to  determine  its  environmental 

hazard.  Graphs  are  given  for  the  fields  at  a  height  z  from  a  HED 

buried  0.175  meters  in  dry  soil  and  operated  at  144  MHz.  But  the  values 

given  for  conductivity  (a  ■  4  x  10  ^  mhos/meter)  and  permittivity 

8 

(eg  ■  4e0)  do  not  satisfy  the  condition  that 


t2-  <  0.05  (ie  n  >  20) 
k«“  *“ 


(2-6) 


where 


ka  -  oo(yoeo)' 


(2-7) 


k  -  (u>2y0e  +  j«MoO‘ 

o  o  o 


(2-8) 


For  these  values  of  a  and  e  ,  the  imaginary  part  of  k  is  at  least  an 

8  8  8 

order  fo  magnitude  smaller  than  the  real  part  for  frequencies  above 
2  MHz.  So  at  144  MHz, 


k  a  (4u0eo)' 


(2-9) 


0.5 


(2-10) 


Therefore, 

k  . 
a 

k  " 

g 

which  is  a  factor  of  10  too  big. 

A  different  approach  is  given  by  Bannister  (Ref  11)  that 
combines  the  Sommerfeld  method  with  image  theory.  He  gives  equations 
for  the  far-field  and  requires  that  n*>15.  He  also  requires  that  the 
observation  distance  be  at  least  three  times  greater  than  the  burial 
depth.  A  few  results  are  given  for  the  HF  range,  but  there  seems  to 
be  no  other  restrictions  on  frequency. 

In  all  the  reports  mentioned  so  far  (Ref  2-11),  the  current 
distribution  on  an  actual  antenna  is  not  analyzed.  Calculations  are 
performed  only  for  an  infinitesimal  dipole, with  a  uniform  current 
distribution.  Entzminger  et  al  (Ref  6.7)  do  replace  the  dipole  moment, 
Idl,  by  a  general  dipole  moment, 


L/2 

Dm  -  I  I(x)  exp(jk  xcos<J>)dx 
-L/2  0 


(2-11) 


but  does  not  evaluate  it.  Also,  this  dipole  moment  is  only  valid  for 
the  fields  at  the  surface. 


Moment  Method  and  Others 

The  moment  method  has  been  applied  to  the  buried  antenna  problem 
to  find  the  pattern  without  knowing  the  current  distribution.  Miller 
and  Burke  (Ref  12)  and  Miller  and  Deadrick  (Ref  13)  use  the  moment 
method  on  the  thin  wire  integral  equations  along  with  a  transmission 
coefficient  for  a  buried  HED.  They  report  their  results  to  be  consis- 
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tent  with  the  Sommerfeld  method.  It  seems  to  be  applicable  to  VHF  but 
the  mathematics  are  outside  the  scope  of  this  thesis.  Numerical 
results  are  given  in  Reference  13  for  a  buried  wire  antenna  operated 

J 

at  frequencies  of  1  MHz  and  10  MHz. 

Brammer  (Ref  14)  presents  the  moment  method  and  a  transmission 
line  method  as  two  ways  to  calculate  the  fields  for  a  buried  insulated 
wire  antenna  operated  at  2  to  16  MHz.  The  methods  and  assumptions  are 
not  described  very  well  but  he  reports  them  to  be  consistent. 

Many  other  references  to  buried  antennas  were  found  but  only  the 
ones  that  are  discussed  here  apply  directly  to  the  problem  of  finding 
the  far-field  pattern  of  a  buried  HED.  Some  of  these  other  references 
are  listed  (Ref  18-43)  for  future  reference. 
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Perfect  Dielectric  Method 


If  the  ground  is  assumed  to  be  a  perfect  or  slightly  lossy 
dielectric  (e'*'«  s'*),  standard  ray-optics  and  antenna  theory  can  be 
used  to  calculate  the  far-field  antenna  pattern  for  a  buried  antenna. 
Table  III-]  shows  the  complex  index  of  refraction. 


ler  Joje0 


(3-1) 


for  different  soils  with  to  -  2ir  x  37.5  MHz.  For  dry  conditions,  the 

imaginary  part  of  n  can  be  neglected  and  the  slightly  lossy  dielectric 

S 

assumption  holds.  For  these  values  and  at  this  frequency,  the  assump¬ 
tion  is  not  very  good  for  wet  soils. 

Assuming  dry  soil,  the  field  distribution  of  the  buried  antenna 
can  now  be  calculated  by  standard  techniques.  The  method  in  Stutzman 
and  Thiele  (Ref  15.25)  is  modified  to  account  for  the  air-dielectric 
interface.  The  E-field  at  some  far-field  point  Q  for  a  dipole  of 
length  L  on  the  x-axis  as  shown  in  figure  III-l  is 


Eg  “  "  4irR 


(cos0  cos<f)0  -  sin$$)Dm  exp(-jk  R)  (3-2) 

8  8 


where 


u 

-  ^  I(x)  exp( jk^xsin0^cos4>)dx 


(3-3) 


k  ■  n  k 
g  go 


(3-4) 


When  the  half-space  above  z«h  is  replaced  by  air,  the  radiation 


observed  at  point  Q  is  different  due  to  the  following: 


Figure  III-l.  Orientation  of  the  HED  in  an  Infinite  Dielectric  Medium. 
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1)  Only  part  of  the  radiation  incident  on  the  interface  is 


transmitted;  the  rest  is  reflected. 

2)  Refraction  changes 

a)  the  direction  of  the  radiation, 

b)  the  dipole  moment,  Dm,  and 

c)  the  radiation  intensity,  U^. 

The  Fresnel  transmission  coefficient  (Ref  16.74),  t,'  is  used  to 
calculate  the  part  of  the  field  that  is  transmitted  such  that 


E 

E 


et  "  t0Eeg 

(3-5) 

m  t  E 

4>t  4>  4>g 

(3-6) 

For  the  0-polarized  E-field, 


2n  cos0 

t  - - S - 8 - 

0  n  cos0  +  n  cos0 
a  a  a  g 

and  for  the  ^-polarized  E-field, 


(3-7) 


2n  cos0 

r  m  ^  8 

<b  n  cos9  +  n  cos0 
*  a  g  g 


(3-8) 


The  effect  of  refraction  is  described  by  Snell's  law 

n  sin0  -  n  sin0  (3-9) 

a  a  g  g 

Using  equations  (3-4)  and  (3-9)  in  (3-3)  results  in 

L 

Dm  ■  f  I(x)  exp( jkoxsin0acos<j))dx  (3-10) 

which  is  the  same  dipole  moment  as  for  an  HED  in  free  space.  The 
ground  has  no  effect  because  all  rays  at  a  certain  angle  travel  the 
same  distance  through  the  ground  as  shown  in  figure  III-2.  Because 
array  theory  is  derived  from  the  free  space  dipole  moment,  this  also 
means  that  conventional  array  theory  can  be  used  in  analyzing  buried 
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antenna  arrays. 


The  phase  factor,  exp(-jk  R)  in  equation  (3-2)  changes  slightly 

8 


y  *  exp 


/ 

l  COS0 


-  jk  (R  -  htan9  s 
o  g 


in0a) 


(3-11) 


The  wave  travels  a  distance  h/cos0  through  the  ground,  but  the  distance 

8 

it  travels  through  the  air  is  shortened  by  htan0^sin0a  as  shown  in 
Figure  III-2.  Equation  (3-11)  is  then  simplified  to 


y  »  exp(-jk  R  -  jk  n  hcos@  ) 
'  o  o  g  g 


(3-12) 


by  using  Snell's  law,  a  trigonometric  identity,  and  equation  (3-4). 

If  the  losses  and  near-fields  are  ignored,  then  the  power  just 


above  the  surface  in  the  solid  angle  dft  must  be  equivalent  to  the 

8 

power  in  the  far-field  in  the  solid  angle  dft  ,  such  that 

cl 


u  (0  ,<|>)dft  U  (0  ,<{>)dft 
a  a  a  tg  e 


(3-13) 


where 


dft  -  sin0  d0  dd> 
a  a  a 


(3-14) 


dfl  ■  sin0  d0  dd> 
g  g  g 


(3-15) 


Differentiating  Snell's  law  yields 


n  cos0 

d0  -  -2 - ^  d@ 

g  n  cos0  a 
g  g 


(3-16) 


and  substituting  this  and  Snell's  law  into  equation  (3-15)  results  in 


n  n  cos0  n^cosS 

dft  -  —  sin0  — - -  d0  d<j>  =  -§ - -  dfl 

g  n  an  cos©  a  y  n^cos0  a 

6  g  g  g  g  g 


(3-17) 


From  equations  (3-13)  and  (3-17) 


n2cos0 

o,(9a.*)  -  ve ,.♦> 

g  g 


(3-18) 
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(3-19) 


From  Poynting  s  theorem,  the  radiation  intensity  is 

—  ft  A 

U  «  Re(^E  x  H  )R2R 


In  free  space  or  air,  E  and  H  are  perpendicular  and  H  =  E/r)o»  therefore 


W*’  ^.hf  r! 


(3-20) 


Thus,  the  far-field  pattern  is 

2  ,  ■  n2cos0 

U  (0  ,4>)  =  E  f  - —  R2  ( 

a  a  Y  2rio  I  gl  n  cos0 

8  8 

So  from  equations  (3-2),  (3-5),  (3-7),  and  (3-10),  the  far-field 


(3-21) 


pattern  for  the  8-polarized  E-field  is 


U  (9 .,♦> 


2n  cos0 


2/  1  \/-jwy 


1  COS0 
g  a 


!f!s _ A7— ¥• 

+  nacos8^\2no/V 


jtuy  \2n2cos0 

cos0  cost})]  — | - 

4tt  g  /  nzcos0 

'  -  g  g 


(i «» 


\2 

)  exp( jkosin0acost}))dx  1 


(3-22) 


where  y  =  y  . 
8  ° 


Results  using  this  equation  will  be  shown  in  the  next  chapter 


after  a  brief  look  at  the  results  from  the  Sommerfeld  method 


IV.  Results 


As  seen  from  the  preceeding  twOi  chapters,  there  are  several  ways 
of  analyzing  the  buried  antenna.  The  next  step  is  to  apply  these 
methods  for  this  specific  case  and  compare  the  results.  It.  would  be 
even  better  to  compare  them  to  actual  measurements,  but  these  are 
unavailable  at  this  time. 

Of  the  three  main  methods;  the  moment  method  (Ref  12,13),  the 
Sommerfeld  method  (Ref  2-11),  and  the  dielectric  method,  only  the 
Sommerfeld  method  and  dielectric  methods  are  discussed  here.  The 
moment  method  looks  promising  but  is  better  left  for  another  study. 

Biggs  and  Swarm  (Ref  4)  and  Vaziri  (Ref  8)  use  the  Sommerfeld 
method  to  give  solutions  that  seem  the  most  applicable  here.  Many  of 
the  other  solutions  were  not  usable  because  of  the  specific  assumptions 
or  conditions.  Usually,  the  analyses  applied  only  for  lower  frequencies 
and  hence,  the  index  of  refraction  is  assumed  to  be  greater  than  what 
it  is  fit  VHF . 

Vaziri 

The  Fortrai)  program  written  by  Vaziri  looked  very  promising  as 
mentioned  before.  It  was  typed  into  AFIT’s  Vax  11/780  computer  line 
by  line  as  given  in  his  dissertation.  Some  of  it  was  excluded  because 
the  code  given  for  the  vertical  electric  dipole  and  magnetic  dipoles 
was  not  needed.  At  first,  the  program  would  not  compile.  No  syntax  or 
error  message  was  generated  during  assembly.  Only  when  the  multiple 
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-V- 


■  v\> 


RETURN'S  were  taken  out  of  the  FNC  function  did  the  program  compile 
without  generating  a  compiler  error.  Evidently,  there  were  more 


RETURN'S  in  the  FNC  function  than  the  Vax  could  handle.  When  these 

i 

were  reduced  by  deleting  more  unneeded  code,  the  program  compiled. 

This  version  is  shown  in  Appendix  A  for  future  reference. 

When  it  was  run  with  n  =  3.0,  a  *  0.005  mhos/meter,  and  the 

S  S 

burial  depth,  h  -  0.5  meters,  the  output  was  obviously  wrong.  For  a 


distance  of  50  meters,  the  antenna  pattern  increases  toward  infinitity 
at  an  angle  of  35°  off  the  horizon.  For  greater  distances,  the  angle 
where  the  pattern  goes  to  infinity  is  decreased.  The  code  was  examined 
several  times  but  the  error  was  not  found. 


Biggs  and  Swarm 

Biggs  and  Swarm  give  the  closed  form  solution  of  the  far-field 
0-polarized  E-field  as 


j  cos0  Jn2  -  sin20 

E,  -  -  60k  Idl  cos*  r - 

bs  R  o  ¥  ^  q  ft I  .  2q 

n  coso  +  4/n  -  sin  0 

g  a  If  g  i 

•  exp(jk  R  +  jk  hJn2  -  sin20  ) 
o  o  ™  g  a 


(4-1) 


They  have  assumed  that  k  R/n  »  1  .  This  is  met  when  R  >  100  meters 

o  g 

for  n  ■  3.0  .  Using  the  same  argument  as  for  the  dielectric  method, 

% 

the  imaginary  part  of  n  can  be  neglected  at  this  frequency.  The 

g 

radiation  intensity  can  be  calculated  from 
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Ubs  ’  HPbs1  R 


2a2 


(4-2) 


(see  equation  (3-20)).  Then  from  Snell's  law  and  a  trigonometric 
identity 


it  can  easily  be  shown  that 


cos@ 


U. 


bs 


U 


a  cos0 


g 


(4-4) 


where  U  is  the  result  from  the  dielectric  method,  equation  (3-22),  with 

a 

dipole  moment  Idl. 

This  seems  to  indicate  that  these  two  methods  of  calculating  the 

G-fields  are  compatible.  The  results  are  in  complete  agreement  at 

0  “0°  (directly  overhead),  but  diverge  as  the  angle  increases  toward 

the  horizon.  For  n  >  3  ,  cos0  -  1  for  all  0  ,  so  the  difference  is 

g  ~  g  a 

mainly  a  factor  of  cos0  .  The  equations  given  in  the  other  reports  by 

cl 

Biggs  and  Swarm  (Ref  2,3,5)  as  well  as  King  and  Shen  (Ref  10)  also 

give  the  same  answer  for  0  *  0° . 
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Computer  Program  for  the  Dielectric  ■Method 

The  Fortran  program  shown  in  Appendix  B  was  developed  from  the 
equations  of  Chapter  III  and  used  to  calculate  the  radiation  pattern  i 
for  a  buried  half-wave  HED.  Although  the  program  will  calculate  both 
0  and  <|)-polarizations,  only  results  for  the  0-polarization  are  shown 
here  for  simplicity.  No  new  information  would  be  gained  by  looking  at 
the  ^-polarization  but  it  might  be  important  in  a  later  study  so  the 
equations  were  left  intact. 

The  integration  to  calculate  the  dipole  moment  (equation  (3-10)) 
is  performed  numerically  because  the  effect  of  burial  on  the  current 
distribution  is  not  known.  In  this  way,  different  current  distributions 
can  be  tested  to  see  what  pattern  is  generated. 

The  Uff  subroutine  (page  B-3)  can  be  used  to  calculate  the  far- 


field  radiation  intensity  in  any  direction  in  the  air  above  a  buried 


,  -  .■%  .'s  4*s  '  , 


HED.  Both  U  and  U.  are  calculated  given  the  index  of  refraction,  0  , 
a  os  a 

and  the  parameters  for  the  current  distribution.  The  program  should 
be  applicable  to  other  frequencies  by  changing  the  value  of  FQ  (frer 
quency)  as  long  as  the  imaginary  part  of  the  index  of  refraction  is 
negligible  at  that  frequency. 


i 


Calulated  Patterns 

The  results  shown  here  are  for  the  radiation  intensity  in  the 
x-z  plane.  A  sinusoidal  current  distribution 

I(x)  -  sin(krkQx)  exp(-ax)  (4-5) 

was  used  to  model  the  current  on  a  4-meter-long  wire  antenna. 

Figure  IV- 1  shows  the  pattern  for  the  half-wave  HED,  both  in 
free-space  and  buried.  The  buried  antenna  radiates  a  much  more  constant 
pattern.  Figure  IV-2  shows  the  change  in  the  pattern  for  the  range  of 
soil  dielectric  constants.  The  radiation  intensity  near  the  horizon 
does  not  change  very  much,  but  these  do  not  take  into  account  how  the 
current  distribution  might  change.  Here,  the  current  distribution  was 
held  constant,  but  in  the  real  world  the  current  distribution  may  change 
drastically. 

Figures  IV-3  -  and  IV-4  show  what  happens  as  the  current 
distribution  is  changed.  In  Figure  IV-3,  the  current  has  been  atten¬ 
uated  by  changing  the  value  of  a  in  equation  (4-5).  This  does  not 
change  the  pattern  except  to  lower  it  as  the  attenuation  is  increased. 

Aa  the  propagation  constant  is  changed,  the  pattern  changes  much  more, 
although  it  is  still  the  same  near  the  horizon  as  shown  in  Figure  IV-4. 
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Figure  IV-2.  Calculated  Pattern  for  the  Buried  HED  With  Different  Soil 


Dielectric  Constants. 
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Figure  IV-3.  Calculated  Pattern  for  the  Buried  HED  With  Different 


Antenna  Attenuation  Coefficients  (e  =  9.0) 
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Figure  IV-4.  Calculated  Pattern  for  the  Buried  HED  With  Different 


Antenna  Propagation  Constants  (e 


9.0). 


Direct ivity  and  Gain 

It  is  more  desirable  to  plot  the  directivity  or  gain  rather 
than  just  the  pattern.  That  would  mean  that  the  whole  pattern  would  be 
shifted  up  or  down  according  to  the  total  amount  of  power  radiated  or 
used.  The  directivity  is  given  by 

D(e,<t>)  -  Hiii&i  (4-6) 

ave 

where 

u  -  —  p  ..  .  -  -r-  u(e,<t>)  dft  (4-7) 

ave  4u  radiated  4it 

To  find  the  total  radiated  power,  the  whole  pattern  must  be  known. 

That  includes  the  radiation  intensity  in  the  ground  that  results  not 
only  from  direct  radiation  from  the  antenna,  but  also  reflections  from 
the  interface.  Then  the  whole  pattern  would  have  to  be  integrated 
numerically.  An  analytic  integration  would  be  even  more  complex  -  if 
it  could  be  done  at  all. 

The  gain  is  given  by 

0(0,$)  .  (4-8) 

P. 

4n  in 

where  is  the  total  power  input  at  the  antenna  terminals.  Both  the 
antenna  impedance  and  the  current  must  be  known  to  find  this.  A  current 
has  been  assumed  but  the  impedance  would  depend  heavily  on  the  type  of 
cable  used,  the  parameters  of  the  ground,  and  the  depth  of  burial. 

As  a  result  of  these  problems,  the  patterns  shown  are  "floating". 


That  is,  their  absolute  values  cannot  be  compared,  but  the  relative 
values  of  one  pattern  can  be.  Therefore,  the  numbers  given  for  the 
magnitudes  of  the  patterns  show  only  the  relative  values  of  a  pattern. 


Losses 


Because  the  soil  is  not  a  perfect  dielectric,  it  will  absorb  a 
portion  of  the  power  as  the  wave  propagates  through.  The  attenuation 
constant  for  a  slightly  lossy  dielectric  is  given  in  Ramo  £t^  al_  (Ref 
17:335)  as 

k 

a  »  v  (4-9) 

where  t'  *  e  £  and  e^=  a  /w  .  The  power  loss  will  then  be 
r  o  g 

Yl  -  (exp(-ad))2  -  Lp(- _|_S_ -j-g-jV  (4-10 

'  go  g  / 

where  h/cos0  is  the  distance  the  wave  travels  through  the  ground.  Bi 
g 

for  n  >  3.0  ,  cos0  -  1  .  Then,  using  k  =.£Tk  ,  and  n  “  k  /we  , 

g“  g  g  \  g  o  *  o  o'  o  ’ 

equation  (4-10)  simplifies  to 


V  h 

Y-.  -  exp( - s=~) 

V  8 


(4-1 1: 


Table  IV- 1  shows  the  loss  for  the  different  soil  parameters  found  in 


the  literature. 


Biggs  and  Swarm 
(Ref  3:206) 

L? 

10 

0.0005 

Entzminger  et  al 
(Ref  6:45-48) 

H 

10  MHz 

60 

0.03 

King  and  Shen 
(Ref  10:1052) 

dry 

144  MHz 

4 

0.00004 

Corvus 
(Ref  :33) 

dry 

HF,  VHF 

9 

0.001 

wet 

HF,  VHF 

25 

0.01 

Burke  et  al 

(Ref  :997 

dry 

HF,  VHF 

9 

0.005 

wet 

HF,  VHF 

20 

0.05 

Table  IV- 1.  Soil  Losses  at  a  Burial  Depth  of  0.5  Meters. 
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V.  Conclusions  and  Recommendations 


Conclusions 

As  stated  in  the  introduction,  a  goal  of  this  thesis  was  to 
develop  a  simple  but  accurate  method  to  determine  the  far-field  pattern 
of  a  buried  HED.  A  relatively  simple  method  has  been  developed,  but 
its  accuracy  has  yet  to  be  determined,  especially  since  a  comparison 
with  measurements  would  be  required. 

The  literature  search  revealed  a  number  of  articles  on  the  topic. 
Most  includes  assumptions  so  that  they  could  not  be  applied  to  this 
case  of  a  HED  operated  at  37.5  MHz. 

The  solution  developed  in  Chapter  III  involves  conventional 

antenna  theory  and  ray-optics.  It  assumes  the  ground  to  be  a  slightly 

lossy  dielectric,  which  was  shown  to  be  a  fairly  good  approximation, 

depending  on  the  moisture  content  of  the  soil.  The  difference  between 

this  solution  and  the  most  applicable  solution  from  the  literature  is 

a  factor  of  cos0  where  0  is  the  angle  off  the  normal  to  the  ground, 
a  a 

So  at  least  overhead,  the  two  solutions  are  in  agreement.  Hear  the 
horizon,  the  difference  is  fairly  large. 

The  directivity  or  gain  was  not  evaluated  because  of  the 
difficulty  of  finding  the  total  radiated  or  input  power.  Several 
patterns  were  shown  for  the  expected  current  distributions.  Burial 
seems  to  lower  the* pattern  overhead  so  most  of  the  power  is  radiated 
close  to  the  horizon.  Most  importantly,  it  was  shown  that  the  pattern 
is  dependant  on  the  dielectric  constant  of  the  soil  and  how  it  affects 


the  current  distribution. 


Recommendat ions 

In  order  to  more  fully  compare  what  happens  when  an  antenna  is 
buried,  the  program  in  Appendix  B  should  be  expanded  to  integrate  the 
pattern  over  the  whole  sphere  surrounding  the  antenna.  This  would  give 
a  value  for  the  total  power  radiated,  and  hence,  the  directivity  could 
be  determined.  Also  it  would  be  interesting  to  find  out  exactly  what 
happens  to  the  current  distribution  and  impedance  on  an  antenna  as  it 
is  buried. 

The  difference  between  the  Biggs  and  Swarm  equation  and  the 
dielectric  method  should  be  examined  to  find  why  they  are  the  same 
overhead  but  different  at  the  horizon.  Finally,  it  would  be  good  to 
get  Vaziri's  program  to  work  and  compare  those  results.  These  tools 
could  then  be  used  to  more  effectively  evaluate  the  buried  antenna  array 
of  Part  II. 
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Appendix  A 


Vaziri's  Program 


jfs  ^  ^  ^  4p  *  ^  ^  *p  ^p  *  'p  *  »p  A*  *  *p  ^  ^  ^  ^  ^  T  ^n  *p  ^  ^  *p  AA  *p  *p  ^  ^  »p  »o  *p  A*  ^  ^  *r  ^  'p  ^  AA  v  *p  ^  ^  A*  *  *  *p  A*  *  t  ™ 

*  Main  program  to  calculate  the  far-field  pattern  # 

*  using  F.  Vaziri's  program  # 

4c  1)  Reads  data  from  'datath'  * 

*  2)  Calls  Vaziri's  solution  for  each  point  in  the  * 

*  desired  pattern  4 

*  3)  Formats  data  for  plotting  # 

###*#*###**#*#****r*#r*###:f<*##*#****#**#********#******* 

real  N,X<5,0J90) , Y<5,0J90> 
data  <  <Y<i,.j>  ,J=0,90)  ,i  =  l , 5)/455*-50 ./ 
open  < 1 ,f ile=' datath ' ) 
rewind  1 

read<l rl)  Nplots, Id 
t  format<2i3> 

30  write<£,2)  Id,  Nplots 

2  format</////'*',i2,5x,i3,'  plots') 
write<6,3) 

3  format </'  pi ot ' , 3x, 'N' ,?x, 'Sg ' ,8x, 'p  ') 

^1^«b  ^  A  «A»  lb  ^  ^  ^  ib  ^  ib  lb  lb  ^  ib  ^  ^  W  ♦  ^  ^  U/  ^  ^  ^  ib  ^  lb  lb  O,  ^  U,  ^  lb  J/  ^  ^  ^  ^  \b  ^  ^  ^  \J/  ^ 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ip  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ip  in  ^  ^  ^  ^  ip  ^  ^  ^  ^  ^  ^  ^  ^  A*  ^  * 

*  do  loop  to  plot  different  cases  on  the  same  graph  # 

jfc  jfc  jfC  ^  jjj  ij^  j|^  ^  ij^  ij^  ^  t|^  ^  ^  «j^  ^  jj^  i|^  ||^  j|^  ^  ^  ^  i|f  ^  ^  ^  ^  ^  ^  ^  ^  ^  i|^  ^  ^  ^  ^  i^i  ^ 


do  10  Ip=l,  Nplots 

read (1,4)  N,p,Sg 
format<3fl0.5) 


write(6,5)  Ip»N,Sg,p 
5  format<i2,5f J0.3,2il0) 

JfC  JfC  J^C  JfC  J^C  *  j|(  )j(  ^  ^  ^  ^  ^  ij^  ^  ^  j|^  ^  j|g  ||^  i|^  ^  i^  i^i  ^  ^  ^  ^  ^  \|r  i|^  ^  ij^  i|^  ^ 

$  do  loop  to  vary  XI  from  0-90  degrees  # 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ip  *  ip  4A  *  ^  ip  ip  ip  ^  ^  A™  *  aA  ip  ip  ip  ^  ip  ^  A*  *  ^  ^  ip  ^  ^  ^ 

SO  do  20  1X1=5,85,5 

Xld-IXI+.OOOl 

call  Pviz<XId ,N,Sg ,Pw,p ) 


Y( Ip , IXI )=Pw 
X< Ip , JXI )=XId 


o.v.nVAc.o.;.  •-•  -•  •  .• 


•*"  *.  \  *■- 


if (Pw»ge»te3)goto  10 


continue 

continue 


write<6,8> 

forir.ot<//J5>i,  '1 '  ,9x,  '2'  ,9x,  '3'  ,9x,  'A'  ,9x,  '5'  ) 
write<£,9)  <X<  1 , 1 ),  <Y(.j  ,  1  >  ,J  =  1 ,5) ,  1=0,90,5) 
format<6f 10*3) 


reQd(l,l)  Nplots,Id 
if <ld,ne»0)goto  30 


stop 

end 


I'llWXMVWW 


zv 


r; 


subroutine  Pviz(XId  »N, Sg  ,Pw,p  ) 

'V  'i  *V'  ^  ^  |1‘ <V ^  *1^  4/  ^  ^  <t  W  ^  ^  W  ^ ^  ^  ^  ^  4  ^  ^  W  ^  'V  ^ ^  ^  'V  ^ 4* ^ 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  'J*  ^  ^  ^  ^  ^  ^r  ^  ®  ^  o*  i|»  a*  ^  ^  ^r  ^  ^  ^  *  *  <p  ^  ^  ^  ^  ip  ^  ^  ^  ^  ip  ^  ^  ^  ip  ^  ^  ^  ip  ^  ^  ^  ^  ^ 

*  this  subroutine  calculates  the  far-field  radiation  intensity  * 

*  using  a  modified  program  from  a  dissertation  by  F*  Oaziri*  * 

«fe  ^u  ^  ^  ^  U>  ^  J/  ^  ^  ^  W  ^  W  V/  \V  ^  U«  tb  ^  ^  ib  ^  d>  lb  U/  ^  \b  V/  ib  \b  ^  ib  U/  ^  ^  ib  lb  ■b  tb  ib  ib  b  lb  b  \V  ^  )l/  ^  ^  \b 

^  ^  ^  ^  ^  ^  ^  »p  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  <p  ip  ^  ip  ^  ^  ^  ^  ^  n*  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ifi  ^  ip  ^  ip  ^  ^  ^  ip  ^  ^  ^  ^  S|V  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 


d<  <b  b  |b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  il,-  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ip  ip  ^  ^  ^  ^  ^  ip  ^  ^  ^  ^  ^  ip  ^  *y* ^  ^  ^  ^  ^  1^*  ^  ^  ^  ^  ^  ^ 


* 

* 

* 

* 

* 

* 

* 


Input  from  main  program 

Xld  angle  off  horizon  (degrees) 

index  of  refraction  for  ground 
conductivity  of  the  ground 
distance  from  antenna 


N 

Sg 

P 

Output 

Pw 


* 

* 

* 

* 

* 

* 

* 


radiation  intensity 

b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b 

^  ^  ^  ^  ^  ^  "I*  *p  Ip  Ip  ^  Ip  Ip  ^  ^  Ip  ^  Ip  ^  Ip  ^  Ip  Ip  ^  1^  ^  Ip  ^  Ip  Ip  Ip  Ip  ^  Ip  ^  ^  Ip  ^  Ip  ^  Ip  ^ 

.  real  N 

DIMENSION  INDEX (6) 

COMMON  NN , NX , XNEAR , X J  NTER , XFAR 

COMMON  /BLK1/H,Z,FG,ER1,STG1,MU1,ER2,SIG2,MU2,SIGN,R,W2MU0,E0I 
common  /b4/Eint 


^p  ^p  ip  ip  ip  ip  ip  ^p  ^p  ^p  ^p  ^p  ^p  ^p  ^p  ^p  ^p 

*  INDEX  indicates  which  polarisation  * 

*  H  burial  depth  * 

*  FQ  operating  frequency  # 

*  ^  ^  ^  ^  Ip  *  ^  Ip  *  *  *  ^  *  Ip  ^  *  "  4*  l|*  *  *  ^  ^  ^  ^  !p  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

data  PI/3* 1415926536/ 

DATA  INDEX/O, 0, 1,0, 0,0/ 

DATA  H*FG/ *175, 144  *E6/ 


t****###********^*##***#****##*!*:#*#**********#******##**)*:** 

*  ER1 ,  SIG1  dielectric  constant  and  conductivity  of  the  * 

*  half-space  that  includes  the  antenna  * 

*  ER2,  SIG2  dielectric  constant  and  conductivity  of  the  * 

*  half-space  opposite  the  antenna  # 

t**#*#*##****#*******#***#**##*###*#**#*#**#****#*#*******# 

ER1=N*N 

SJGl=Sq 

DATA  ER2,SIG2/1 *0,0*0/ 


CALL  CONST 


Z*#*#**#*****#**#**#*# 
*  convert  deg  to  rad  # 
*********** *********** 
XIr»XId*PI/180. 


***************************************** 


*  Z 

*  R 

* 

* 


vertical  observation  distance 
horizontal  observation  distance 
the  cordinate  system  is  upside- 
down  because  H  must  be  positive 


* 

* 

* 

* 


A-3 


i|f  ^  ^  ^  ^  ^  ^  U/  ^  ^  dr  ^  U/  ^  ib  ib  ^  ^  ib  ^  ib  ib  ^  di  ^  |b  <b  ^  |b  b  ib  ^  « i>  |b  |b  ib  d* 

^  ^  ^  ^  ^  ^  ^  V  V  ^  ^  ^  ^  T"  "P  *  <P  *  *  ^  ^  ^  ^  ^p  *  *  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

Z=-p*sin (XI r ) 

R=p#cos<XIr) 

i 

j  DO  4  NN=1 f6 

A  !  CALL  COMP < INDEX <NN>> 

ib  ^  ^  Ui  ^  dr  ^  ^  ^  ^  ^  ^  4  W  ^  ^  vb  ib  u>  d' 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^p  ^  ^  ^  ^  *p  »p  *  ^  *  »p  ^  ^  /p  rp  ^  ^ 

*  calculate  radition  intensity  * 
Pw-Eint*Eint/< 2*377) 


A 


return 


.Wi 


1 


SUBROUTINE  CONST 

THIS  SUBROUTINE  CALCULATES  NECESSARY  CONSTANTS. 

• 

COMMON  /BLK1/H , Z ,  FQ ,  ER1 , SIG1 , Mill  ,  ER2 , SIG2 , MIJ2 , SIGN , R ,  W2MIJ0 , EO 
COMMON  /BLK2/W, KO, NISO, N2SQ, El, E2,K1,K2,K1S,K2S,N,CSIG, POLE, P 
C  MUO,EG, ABSN 

COMMON  /BLK4/AI,AIK1,AIK2 
COMMON  /BLK6/  KR,E1INV 

COMPLEX  N1SQ,N2SQ,E1,E2,NJ ,N2,K.l  ,K2,N,CSIG  ,KJS,K2S 

COMPLEX  KR,E1INV 

REAL  K0,MIJ0,MU1,MU2 

E0=8.854E-12 
MUO=l .257E-06 
PI=3. 14159 

i 

MU1  AND  MU2  ARE  RELATIVE  PERMEABILITIES  OF  MEDIUM  1  AND  2 

MU1«1. 

MU 2=1 . 


IF(2)  1,2,2 
SIGN=-1 
GO  TO  3 
SI6N=1 ♦ 
CONTINUE 


W=2.*PI*FQ 

C«1./<E0*MU0>**.5 

W2MU0=W*W#MU0 

EOINV*l./EO 

KO=*W/C 


N1SQ*CMPLX<ER1 , SIG1/W/EO) 

E1«E0#N3SQ 

N1*CSQRT  <N1SQ) 

K3  *KO#Nl 
AIK1=AIMAG<K1> 

K1S-K1*K1 

N2SQ=CMPLX <  ER2 , SI G2/W/E0 ) 

E2»E0*N2S0 

,N2*CSQRT(N2SQ> 

K?=K0*N2 

AIK2=AIMAG<K2) 

K2S*K2*K2 


v.v.v, 


mJi  -  -  ■_  -  -  - 


A- 5 


KR=E1/E2 

N=K2/K1 

ABSN=CABS<N> 

CSIG=CMPLX<SIGJ.,-EO*W> 

P0LE=AIMAG<C5GRT<  <E2*E2*K1S-E1*E1*K2S)/<E2#E2-E1*E1 > ) ) 

RETURN 

ENB 


SUBROUTINE  COMP (INDEX) 


THIS  SUBROUTINE  EVALUATES  THE  FIELD  COMPONENT  IF  INDEX  NE  0 
COMMON  NN , NX , XNFAR f  X I NTER , XFAR 

COMMON  /BLK1/H,Z,FQ,ER1 ,SIG1,MU1 ,ER2,SIG2,MU2,SIGN,R,W2MUP0,Ei 
COMMON  /BLK2/W ,  KO ,  N.t  SO ,  N2SO , El , £2 , K1 , K2  f  K1 S ,  K2S ,  N ,  CSIG  f  POLE  ,  P 
MUOf EO» ABSN 
COMMON  /BLK3/CDNST. 

COMMON  /BLK4/AI  »AIK1»AIK2 
COMMON  /BLK5/IFLG1 ,IFLG2 
common  /b4/Eint 

COMPLEX  N1SG,N2SQ,E1 ,E2,K1 ,K2,N,CSI6rKlSrK2S 
COMPLEX  I 1, 12 , I, CONST 

REAL  KOrMUO 

IF  (INDEX. EQ.O)  GO  TO  5 

IFLG1=0 

IFLG2=0 

EPS=0. 

I1»<0.,0.) 

I 2= (0**0.) 

START=0. 

END=.9?*P0LE 

INTEGRATE  FROM  0  TO  THE  VICINITY  OF  THE  POLE 

CALL  AUTO (START * END* EPS *11*12) 

START=END 
END=1 .005*P0LE 


INTEGRATE  ON  THE  POLE 

CALL  AUTO (START *  END* EPS* I 1*12) 

START=END 

END*1 .0.1 5#P0LE 

IFLG1=0 

IFLG2=0 

CALL  AUTO (START, END* EPS* Ilf  12) 

IFLG1*0 

IFLG2=0 

START=END 


CONTINUE  INTEGRATION 


DO  2  J=1 » 41 » 4 


END=P0LE*(2.**AJ) 


TERMINATE  INTEGRATION  IF  INTEGRAND  IS  NEGLIGIBLE 


IF  (IFLG1+IFLG2-2)  1,3,3 
CAL  L  AUTO ( START , END , EPS ,11,12) 
START=END 
CONTINUE 
1=11+12 

Eint=cabs<  I *C,ONST ) 


RETURN 


SUBROUTINE  AUTO < START , END , EPX , I J ,12) 


THIS  SUBROUTINE  PERFORMS  NUMERICAL  INTEGRATION  ALONG  UERTIC 
BRANCH  CUTS* 

COMMON  NNrNX 

COMMON  /BLK1/H » Z  *  FQ » ER1 >  S JG1 » MU1 > ER2 * SIG2 f MU2 f  SIGN W2MU0 y EO 
COMMON  /BLK4/AI,AIK1 „AIK2 
C0MM0N/BLK5/JFLG1 , IFLG2 

COMPLEX  I1,I2,T1,T2,CSIQD,FNC 

EXTERNAL  FNC 

ETA  IS  INTEGRATION  ERROR 

ETA=.01 
NX=  1 
AI=AIK1 

IF<  <2+H) .LT .0* )  NX=3 

TERMINATE  INTEGRATION  ALONG  BRANCH  CUT  1  IF  INTEGRAND  TOO  S 

IF  (IFLG1.EQ.1)  GO  TO  1 
T1=CSIGD  <  START , END , FNC , EPX , ETA) 

IJ=I1+T1 

IF  <CABS<T1 ) »  LT. CABS<I1>*. 01 >  IFLG1=1 
CONTINUE 

NX=NX+1 

AI-ATK2 

TERMINATE  INTEGRATION  ALONG  BRANCH  CUT  2  IF  INTEGRAND  TOO  S 

IF  ( IFLG2 .EQ » 1 )  GO  TO  2 
T2*CSI0D  <  START , END , FNC , EPX  ,  ET A ) 

I2-I2+T2 

IF  (CARS<T2) ♦LT*CABS<I2)#*01)  IFLG?=1 
CONTINUE 

EPX*CABS < 1 1 +12  >  * . 001 
RETURN 


COMPLEX  FUNCTION  FNC<X) 


THIS  FUNCTION  SUBPROGRAM  EVALUATES  THE  NUMERICAL  VALUE  OF  TH 
INTEGRAND.  WHEN  ARGUMENT  IS  ZERO,  IT  ALSO  CALCULATES 
CORRESPONDING  BANOS ' S  APPROXIMATIONS. 

dimension  JFLG<6) 

COMMON  NN , NX, XNEAR , X INTER , XFAR 

COMMON  /BLKJ/H,7,FG,ER.1  ,SJG.t,MUJ  ,  ER2,  S1G2,  MU2 ,  SIGN ,  R ,  W2MU0 ,  EOI 
COMMON  /BL.K2/W , KO , NtSQ , N2SQ , El  ,E2 , K1  , K2 , K1S, K2S , N , CSIG , POLE ,  PI 
MUO , EO , ABSN 
C0MM0N/BLK3/C0NST 
COMMON  /BL.K4/AI  ,A.TK.t ,  AIK2 
COMMON  /BLK6/  KR,E1INV 

COMPLEX  KR , KRA2 , KRB2 , E 1 I NV 

COMPLEX  N1SQ , N2SQ ,E1,E2,K1,K2,N,CSIG,K1S,K2S 
COMPLEX  B 1  ,  B2 ,  L. ,  ARG  ,  EXP 1  ,  EXP2 , FRAC1  ,  FRAC2  y  A1 ,  A2 
COMPLEX  XFRAC1,XFRAC2, CONST 
complex  T 

COMPLEX  HONEO,HONE1 
COMPLEX  HO,H.1  ,H.1D,FX.1  ,FX2 

REAL  KO , MUO , MU1 , MU2 

DATA  JFLG/I ,1,0,1, 1,-1/ 


IF(X.GT.l .E-30)  GO  TO  1 
FNC=<0.,0.) 

CALCULATE  OBSERVER  HEIGHT  FOR  BANOS'S  APPROXIMATIONS. 

Z1=Z-H 

SKIP  BANOS'S  APPROXIMATIONS  IF  SOURCE  DIPOLE  IS  IN  AIR 

4 

IF  (ABSN.GT .1 . )  GO  TO  180 

T»  <  0 ♦ f 1 . ) *N*CSGRT ( . 5*P I *  <  0 . , 1 , > *K2*R ) 

IF <  <Z+H) .LT.O. )  GO  TO  72 

IF  <R#< AI+X) .LE.20. )  GO  TO  (3, 4, 6,7), NX 

FNC=<0» ,0. ) 

RETURN 

EVALUATE  FIELD  INTEGRAL  Al.  ONG  BRANCH  CUT  1  FOR  POINTS  OF 
OBSERVATION  IN  MEDIUM 


Bl=CSQRT<X*X-<0.  ,2.  >*I\3  #X> 


print*,*ot  3* 

B2=CSGRT(K2S-K:tS+X*X-<0.  ,2,  >*K3*X) 

L=Kl+(0»  , 1  * )*X 

C  USE  MODIFIED  FORMULAS  IF  SOURCE  DIF’OLF  IN  AIR 

IF  < ABSN.GT . 1 . )  B2=-B2 
£XP.t=CEXP (  (0 . , -3  *  )#B1#SIGN*7> 

EXP2=CEXP (  <0.  ,  1 .  >*B1*<7:-F2.*H> ) 

KRB2=KR*B2 

ARG=L*R 

IF  (JFLG(NN))  200  ,  2.1  0  ,  220 
220  HO=HONEO < ARG ) 

Hl=HONE.'t  <  ARG ) 

HtD-HO-Hl/ARG 
FRAC1= ( B3 +B2)/ <B1-B2) 

-  FXI=FRAC1*EXP2 
210  FRAC2=  <  B1 +KRB2 ) / < B1-KRB2 ) 

FX2=FRAC2*EXP2 
GO  TO  230 

200  FRACl= ( B t +B2 )/(Bl-B2) 

FX1=FRAC3*EXP2 
230  continue 

C  EVALUATE  FIELD  INTEGRAL  ALONG  BRANCH  CUT  2  FOR  POINTS  OF 

C  OBSERVATION  IN  MEDIUM  1 

4  A3=CSQRT(K1S-K2S+X*X-(0,,2,)*K2*X) 

print*, "at  4* 

A2=CSQRT  <  X*X- <  0  * ,2  * ) *K2*X ) 

L=K2+(0.,1.)*X 

C  USE  MODIFIED  FORMULAS  IF  SOURCE  DIPOLE  IN  AIR 

IF  (ABSN.GT ♦ 1 , )  A1=-A1 
EXP2=CEXP(  <0. ,1. >*A1*(7+2.*H> ) 

KRA2=KR*A2 

AR6SL^R 

IF< JFLG(NN ) )  300,330,320 
320  H0=H0NE0(ARG> 

H3=HONE3 (ARG) 

H1D=H0-.H1/ARG 
FRAC1=(A3  +A2)/( A3 -A2 ) 

310  FRAC2=( A1+KRA2 )/( A1-KRA2 ) 

GO  TO  330 

300  FRAC1=( Al+A2)/( A1-A2) 

330  continue 

C  EVALUATE  FIELD  INTEGRAL  ALONG  BRANCH  CUT  1  FOR  POINTS  OF 

C  OBSERVATION  IN  MEDIUM  2 

A  Bl-CSQRT (X*X-<0. ,2. >*K1*X) 

B2=CSGRT ( K2S-K3 S+X*X- <0»,2»)*K1*X) 


L=K1+(0.,61.)*X 
L=Kl+(0.,3 . >*X 


USE  MODIFIED  FORMULAS  IF  SOURCE  DIPOLE  IN  AIR 

IF  ( ABSN.GT .  1 .  >  B2=-B2 
EXPt=CEXP< (0. ,1. >*Bt*H> 

EXP2=CEXP< (0. ,3 .  )*B2*<7-i-H> ) 

KRB2=KR*B2 

ARG=L*R 

IF  ( JFLG(NN) )  400,410,420 
HO=HONEO  <  ARG ) 

Hl=HQNE.t  <  ARG) 

H1D=H0-H1/ARG 
FRAC1=1  ./<  BIT B2) 

FRAC2=1 ./(B1-B2) 

XFRAC1=E  3 INV/ ( B3  +KRB2 ) 

XFRAC2=:E1  INV/  (  B1 -KRB2 ) 

GO  TO  430 
FRAC1=1  . /(  B1. +B2) 

FRAC2=1 */(Bl-B2) 

GO  TO  70 

> 

EVALUATE  FIELD  INTEGRAL  ALONG  BRANCH  CUT  2  FOR  POINTS  OF 
OBSERVATION  IN  MEDIUM  2 

Al=CSGRT(NlS-K2S+X#X-(0» ,2. >*K2*X> 

A2=CSQRT  <  X*X- <  0 . , 2 . ) *K2*X ) 

L=K2+(0.,1.)*X 

USE  MODIFIED  FORMULAS  IF  SOURCE  DIPOLE  IN  AIR 

IF  (ABSN.GT. 1.)  Al=-A3 
EXPt*=CEXP(  (0. ,1. )*A1*H> 

EXP2=CEXP< (0. ,1 ♦ >*A2*<Z+H> ) 

KRA2=KR*A2 

ARG=L*R 

IF < JFLG(NN) )  500,510,520 
HO=HONEO ( ARG  > 

H1=H0NE3 (ARG) 

HID=H0-H1/ARG 
FRAC1  =  1 ♦/( A3  +A2 ) 

FRAC2=l’»/(  A1-A2) 

XFRAC1=E1 INV/ ( A1+KRA2 ) 

XFRAC2=E1 INV/ ( A1-KRA2 ) 

GO  TO  530 
FRAC1=1 »/ ( A1TA2) 

FRAC2=3 ♦/( A1-A2) 

GO  TO  71' 

*  ^  ^  "r  *  ^  ^  ^  ^  ^  ^  ^  4*  ^  ^  ^  ^  ^  ^  ^  4"  ^  ^  4*  *  ^  ^  Ip  ^  ^  *p  »p  ^  ^  ip  ^  ^  ^  ^  ^  ^ 

HED  Ez2 

##**  #  ##*#  *  *>X  *  #*#  *  #  #  *  *  X  **  *  *  #  ###  :«  Xofc  ## 


2 


FNC=L*L*H0N£l<ARG)*EXP2*Bt*<XFRACl/EXPl-XFRAC2*EXPJ.> 

RETURN 

FNC=L*L*H0NE1 ( ARG  >  *A 1 *EXP1 * < XFRAC2*EXP2-XFRAC 1/EXP2 ) 
RETURN 

CONST- 1 ♦  /  (4  ♦  ifcPIifcW) 

JF( ABSN »GT» 1 »  >  RETURN 

XNEAR  =CABS<K1/(2*#PI#CSTG#R#R) ) 
XINTER=CABS(K2*K2/<2.*PJ*CSIG*R)/N*n  ,+T) ) 

XFAR  =CABS(K1/<2**PI*CSIG*N*N*R*R) > 


RETURN 


COMPLEX  FUNCTION  CSTGIK A,B,FCN, EPS, ETA) 

ROMBERG  INTEGRATION  FOR  COMPLEX  FUNCTION  WITH  REAL 
ARGUMENT,  MODIFIED 


COMPLEX  FCN ,Q(11) ,T > SUM , FCNXI , QX1 , QX2 

EXTERNAL  FCN 

H=(B-A)/2. 

T=H* < FCN < A ) +FCN < EO  ) 

NX=2 

DO  12  N=1 , 10 
SUM=0 » 

DO  2  1*1, NX, 2 

FCNXI =FCN ( A+FLOAT  < I >  *H ) 

SUM=SUM+FCNX I 
T=T/2.+H*SUM 
Q<N)  =  (T+H#SUM)/1  .5 
IF<N-2>  10,3,3 
F=4. 

I=N 

DO  4  J=2,N 
1=1-1 
F=F#4 ♦ 

Q<I)=Q<I  +  l)  +  (Q<I  +  l)-a<I) )/(F-l, ) 

IF  <N-3)9,6,6 

X1=ABS  <  REAL  <  Q< 1 ) -0X2) ) +ABS ( REAL ( QX2-QX1 ) ) 

X2=ABS ( AIMAG (0(1) -0X2  > ) +ABS ( AIMAG ( 0X2-QX 1 ) > 

TABS1=ABS ( REAL (0(1))) 

TABS2=ABS (AIMAG (0(1))) 

IF(TABS1 )7 ,8,7 

IF ( XI /TABS 1 -ETA >20,20, 8 

IF(X1-EPS)20,20,7 

TF(TABS2)71 ,81,71 

IF(X?/TABS2-ETA) 11 ,11,81 

TF(X2-EPS) 11,11,9 

QX1=QX2 

QX2=Q( 1 ) 

H=H/2 » 

NX=NX*2. 

URITE(6, 100)  A,B 

FORMAT (43H  ACCURACY  LESS  THAN  SPECIFIED  VALUES — CSIOD  2E11.4) 
CSTQD=Q( 1 ) 

RETURN 


X 


COMPLEX  FUNCTION  HONEO(XZ) 

HANKEL  FUNCTION  OF  FIRST  KIND,  ORDER  0,  IMAGINARY  PART  LT  30 

IMPLICIT  COMPLEX  (X) 

ATXZ=AIMAG<XZ> 

IF (AJXZ »GT *20  ♦  >  GO  TO  83 
XW=-XZ*XZ/4» 

PI=3* 141592653589 
TPI=2  *  *PI 
XI=<0.fl.) 

IF (CABS< XZ >  »ST»5* )  GO  TO  7 
FAC=0. 

XFAC=< 1 ♦ »0» ) 

XSUMY=<0. ,0. ) 

XSUM.J=<  1 .  >0*  ) 

DO  5  1=1 ,50 
AI=I 

FAC=FAC+1 */AI 
XFAC=XFAC*XW/<AI*AI> 

XTE=FAC*XFAC 

XSUM7=XSUMY +XTE 

XSUMJ=XSUM J+XFAC 

IF ( CABS ( XTE >»LT*l»E~tO)  GO  TO  6 

CONTINUE 

WRTTE<6,92)XZ,XTE,XFAC 

FORMAT <5X,'H0NE0  SUM  EXCEEDS  50  TERMS' r6E15»3//) 

H0NE0=<0.f0.) 

RFTURN 

XJ=XSUMJ 

HONEO=X J+XI*  < 2 . /PI ) * ( < CLOG  <  XZ/2 ♦ )  +  ♦ 5772156649 ) *X J-XSUMY ) 

RETURN 

XP=l.-<9./<128.#XZ*X7>)*<3 ♦ -1225 ♦  /< 768 .*XZ#XZ> ) 

XG=-<1*/<8.#XZ> >*(1 .-225./<384.*XZ*XZ) ) 

RZ=REAL<XZ> 

SIGN=1 • 

IF<RZ»LT  *0 ♦ )SIGN=-1 « 

IF<ABS(R2>.LT.TPI)  GO  TO  73 

RZ=R7-SIGN*TPI 

GO  TO  71 

H0NE0=<XP+XI#Xa)#<C0S<RZ-,25*PI>+XI#SIN<RZ+.25#PI>) 

1  *EXP<-AIXZ)*CSC}RT<2,/<PI*XZ>> 

RETURN 

H0NE0=<0*,0.) 

RETURN 

END 


COMPLEX  FUNCTION  HONEl<XZ> 


HANKEL  FUNCTION  OF  FIRST  KIND,  ORDER  1,  IMAGINARY  PART  LT  30 

IMPLICIT  COMPLEX ( X ) 

ATXZ=AIMAG  <  XZ ) 

IF (AIXZ  »GT  *20  » )  GO  TO  81 
XW=-XZ*XZ/4. 

PI-3,141592653589 
TPI=2.*PI 
XI=(0.,1, ) 

IF(CABS<XZ) »GT,5, )  GO  TO  7 
FAC=-. 5772156649 
XFAC=< 1 , ,0, ) 

XSUMY=FAC+ ,5 
XSUM J=  < 1 , , 0 . > 

DO  5  1=1,50 
AI  =  I 

FAC=FAC+1  ,/A.T 
XFAC=XFAC*XW/(AT*<AH-1.  >> 

XTE=  <  F  AC+ . 5/ <  A I  + 1 , ) )*XFAC 

XSUMY=XSUMY+XTE 

XSUMJ=XSUM J+XFAC 

IF<CABS<XTE) , LT, 1 »E-10>  GO  TO  6 

CONTINUE 

WRITE<6,92)XZ,XTE,XFAC 

FORMAT <  5X  »  ' XHFR1  SUM  EXCEEDS  50  TERMS' ,6E1 5,3//) 

H0NE1=<0» ,0. > 

RETURN 

XJ=XSUMJ*XZ/2. 

H0NE1=XI*  <  <  CLOG  < XZ/2 , ) *X  J-l . /XZ ) # < 2 , /PI ) -X7#XSUMY/PI ) +X J 
RETURN 

XP=1.+<15./(128.*X7*X7> )#<1 ,-315,/<256.*X7*XZ>) 

XQ=<3,/<8.*XZ> >*< 1.-35. /<128,*XZ*XZ>> 

RZ=REAL(X7> 

SIGN=1 » 

IF  <RZ.LT ,0, )SIGN=-1 . 

IF  <ABS<RZ) .LT.TPI )  GO  TO  73 

R7=RZ-S.TGN*TPI 

GO  TO  71 

HONE.t  =  (XP+X J#Xfl ) * < COS < RZ- . 75#PI  )+XI#SIN<R7-.75#PI )  )# 
cEXP < -AIXZ ) *CSQRT  <  2 . / <  P I *XZ )  ) 

RETURN 

H0NE1=<0.,0.) 

RETURN 

END 


Appendix  B 

Dielectric  Method  Program 
program  vertic 

****************************************************************** 


*  Moin  program  used  to  calculate  the  far-field  pattern  # 

#  of  a  horizontal  electric  dipole.  # 

#  1)  Reads  data  for  current  distribution  and  soil  parameters  # 

*  2)  Calls  Uff  subroutine  to  calculate  the  radiation  * 

t  intensity  at  each  point  in  the  pattern.  # 

%  3)  Formats  data  for  plotting  # 

$  (plotting  subroutine  not  shown)  # 


#*##*##*###*************  r******************* ********************** 

real  N,Kr,X(5,0t 180)  ,UdB ( 5,0 J 180)  ,U(2) 
common  /b 1/ALP , Kr , shi ft 
-  data  PI/3. 1415926536/ 

open  (l,file='datath'  ) 
rewind  1 

read(l, 1)  Nplots, Id 

1  format(2i3) 

open (2 r file='outth '  ) 

30  write(2,2)  Id,  Nplots  » 

2  format (/////'*' ,i2,5x,i3, '  plots') 

write(2,3) 

3  format(/'  plot'  ,6>c,  'Meth  '  ,6>:,  'J'  ,9k,  'N'  ,9x,  'B'  ,9x,  'ALP'  ,7x 
c,'Kr'  ,8x,'PHId'  ,8x,  'shift'/) 

*  initialize  * 

*  UdB  power  density  pattern  (in  dB's),  * 

*  8err  greatest  integration  error  per  graph  * 

do  50  i=0,180 
do  50  .j=l,5 
UdB(J,i)=0. 

50  continue 

Gerr*le-5 

-a.  WWaA*  a4a  aAa  ^b  ^  ^a  ^  ^  ^  ^  ^  ^  ^  ^  d.  d.  d.  di  ^  ^  ^  ^  ^  ^b  ^  ^  ^  lb  ^  ^  lb  lb  ^  ib  ^  ^  d.  db  ^  ^  ^  lb  ^  lb  ^  ^  ^  ^ 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *  *  *  *  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  * 

*  do  loop  to  plot  different  cases  on  the  same  graph  * 

di  ^b  da  ^b  ^  ib  ^  ib  ^  ^  ^  ib  ib  ^  d.  di  ^  ^  d,  di  d.  d*  >b  ib  d>  <b  yb  di  ^  ^  ^  ^  ^  ib  ^b  d.  ^  lb  ^  ^  ^  da  ^  ^  ^b  ^b  dj  ib  ^b  ib  ib  ^b  ^b  ala 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  )p  ip  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

do  10  Ip  =  l,  Nplots 

4 

read (1,4)  Meth, J,N,B, ALP, Kr,PHId, shift 

4  for»at(2i3,6f 10.5) 


write (2, 5)  Ip,Meth, J,N,B,AlP,Kr,PHId, shift 


5 


format < i5,2il0,6fl0»3) 


yj  ^  ^  «l|*  iL  ^  ^  di  ^  ^  ^  ^  ^  ^  tU  Jfl  ^  y«  dt  dj  ^  ^  ala  lb  ^  ^  y*  ^  ^  ^  ^  lb  ^  ^b  ^  ^  ^  dj  df  lb  lb  da  ^  ab  da  da  ^  ala  ^  ^  ^ 

A  A  A  A  *  ^  ^  ^  ^  ^  »p  ^  H*  ^  ap  ^  ^  ^  ^  ^  ^  *p  T*  ^  ^  ^  »p  ap  ^  Ip  ^  fla  1^  ^  «p  ^a  ^p  ^  ^a  ^  ^  ip  ^  ^  fla  ^ 

*  doi  loop  to  vary  Xld  <the  angle  off  the  horizon)  * 
t  'from  0-90  degrees  for  the  front  lobe  and  then  # 
t  '90-0  degrees  for  the  back  lobe  * 

^da  ^  ^  ^  ab  tb  ^  ^  ib  ill  ib  rb  ab  (b  ib  ^  lb  da  ib  lb  da  ab  ^  ala  ab  *La  b  da  db  ij a  b  lia  tb  via  b  b  b  b  b  b  b  da  b  b  b  b  b  b  b  b 

A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A 

do  20  IXId=0 , 180 
XId=IX.Td  +  *0001 
if  < IXId »gt .90)  then 
PHId=180. 

XId=180-IXId+. 0001 
end  if 

call  Uff <PHId,XId,N, J,B,Err,U< 1 > ,U<2) ) 

-  UdB(Ip,IXId)=10.*alogl0(U<Meth)+le-5) 
X(Ip,IXId)=IXId 

if (Err .gt .Gerr)  0err=Frr 
20  continue 

10  continue 


write(2,ll)  GERR 

11  format(/  greatest  integration  errar=' ,f 10*5///) 

^  j|^  j|^  ^  |||  ||^  ^  y|i  j|^  ^  ^  ^  ^  ^  t ^a  v^a  ^  ^  \|a  i|^  |||  ^  jj^  ||r  i|^  ^  ^  ^  i|^  ^  i|^  ^  j^r  ^  ^  j|^  ^  j|^  ^ 

*  MPL0TS  plotting  subroutine  * 

*  HDC0PY  transfers  plot  to  paper  using  Textronix  4631  # 

*  Hard  Copy  Unit  # 

y  j|^  j|^  ^  ^  ^  a^  ^  a|^  v^a  ^  ^  ^  ^  y^  j|^  ^  ^  ^  ^  ^  y|^  ^  ^  ^  a^  a^  ^  y|^  j|^  y^f  j|^  ^ 

call  MPL0TS(5, 181 , X,UdB,UdBmax) 

write(6,6)  Id 
6  format <i3) 
call  HDC0PY 

yial>  b  •£■  d*  dt  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b 

A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A 

A  test  for  next  set  of  data  * 

^  j|^  j||  j|^  j||  j|^  jj^  j|^  ^  ij^  yj^  ^  y^  y^  yj^  yj^  y^  y^  y^  y^  a^a  y^  y^ 

read < 1,1)  Np lots, Id 
if (Id.ne.O)goto  30 


40 

X 


stop 

end 


B-2 


iL 


subroutine  Uff (PHId ,XId ,Ng , J,  B,  Err,  Upd ,  libs) 

^  ^  «jL»  ^  ^  ^  ^  ^  ■  u  ^  *u  «y  «L  ^  iy  ^  § t>  &y  iy  ^  ^  y  tii  u*  ^  ■!■  ^y  ty  ^  ^  ^  ^y  >y  ty  >y  >y  ^y  <y  iy  ^  ^  «y  «y  ^  » y  y  ^  ^  ^y 

f  f  ^  ^  ^  ^  ^  ^  ^  ^  *  V  *  ^  ^  ^  ^  ^  ^  ^  ^  T  *  *  ^  ^  ^  ^  ^  ^  ^  ^  T  ^  *T  vr  ^  v  ^  ^  ^  *r  *  ^  * 

t  Input  # 

t  B  length  of  antenna  (meters)  * 

*  PHId  azimuthal  angle  (degrees)  # 

*  J  polarization  index  *  * 

$  J=1  for  theta  polarization  * 

$  J=2  for  horizontal  ’  # 

$  Xld  angle  off  horizon  (degrees)  # 

*  Ng  index  of  refraction  for  ground  # 


*  Output 


t  '  Upd  radiation  intensity  (perfect  dielectric  method)  * 

*  Ubs  radiation  intensity  (Biggs  and  Swann  method)  $ 

ft  Err  max  integration  error  per  case  ft 

*  Output  to  INT  * 

*  THafPHIr  * 

SQ(  #  )|(  ^  ^  jj^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  j||  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ||^  ^  |j^  ^ 

external  INT 
common  /b2/THa,PHIr 
real  Ng,Na,E(2> ,Ko 
complex  CDM 

^  ^  ib  ^  ^  ^y  ^  ^  >y  ^  Ut  ^  >y  iy  Ur  ^  «y  <y  ^y  iy  «i.  ty  ^  >y  >y  iy  (y  ^y  ^  ^  ^  iy  <y  ^  <y  ■  y  ^  ^y  ^y  u.  ^y  ^y  ^y 

^  ^  ^  ^  ^  ^  ^  ^  *r  *  *p  ^  ^  ^  ^  ^  ^  ^  ^  ^  *  ■p  ^  ^  *  *p  ^p  *  *  ^  *  ^p  *  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

*  Na  index  of  refraction  in  air  ft 

*  FQ  operating  frequency  ft 

ft  Eo  permittivity  of  free  space  ft 

*  C  speed  of  light  (m/s)  ft 

*  U  radial  frequency  ft 

*  Ko  propagation  coefficient  of  free  space  ft 


*  ETAo  intrinsic  impedance  of  free  space  ft 

•A*  ^  Ai  ^y  <y  ^y  u<  df  «y  a*  iy  iy  ^y  <y  ^  «y  ^y  iy  ^y  ^y  >y  ^  ^y  ^y  u^  *L»  ^y  ^y  x  ^y  ^y  ^  ^y  4y  >y  ^y  ^y  >y  *y  ^y  >y  ^y 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *p  ^  *  *p  ^  ^  *p  ^  ^p  ^  ^p  ^p  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

,  data  Na, PI, FQ,Eo, C/3 .0,3. 1415926536, 37.5eA,8.854e-12,3e8/ 
y«2*PI*FQ 
Ko-W/C 
ETAo=377. 

tttttt*#******x*********#**t 

*  convert  deg  .to  rad  * 

A*  tMm  ib  Ui  a,  y,  a#  dt  ^y  u«  ^y  *y  y«  ^y  ^y  ^y  ^y  ^y  Uj  ui  ^y  »y  ^y  uj  *y  u* 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *  *  *p  ^  ^  ^  ^  ^ 

PHIr=PHId*PI/180. 

XIr«XId*PI/180. 


*****ft***»ft*****ftft#*ft***************** 
ft  snell's  law  * 

*  THa  angle  from  norm  in  air  * 

*  THg  angle  from  norm  in  ground  * 
**»****#**ft**ft»ft»*ftft*ft##ft**ft»ft#ftftft*ft** 

THa*PI/2-XIr 


s 


THg=asin<Na*sin(THa)/Ng> 

i 

*|a  Oi  tb  d<  d>  U#  tb  ^  tb  iji  ^  yb  J»  tb  pb  ^  tb  *1*  ^  tb  ib  ^  W/  ^  ^  tb  pb  tb  b 

^  ^  ^  ^  ^  ^  ^  ^  *T  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

*  transmission  coefficients  * 

*  TT  theta  polarisation  # 

%  TH  horizontal  polarisation  * 
MtXtMMMtXMXXXWWXXMXXXXKX* 

TT=2 . *Ng*cos ( THg ) / <  Ng*c  os ( THa ) +Na*cos ( THg ) ) 
TH=2  ♦  *Ng*cos  ( THg  )  /  <  Ng *cos  ( THg )  +Na*cos  <  THa ) ) 


*#*########*********##***#********* 

#  integrate  to  find  dipole  moment  # 

#  DM  -  dipole  moment  # 

d«  ib  U/  U/  \l/  ^  4'  \l'  v)/  ^  ^  ']/  vi/  vl/  \l/  \1’  'V  ^  ^  <J.*  ^  ^1/  ^  ^ 

J(C  ip  JJC  J|C  3|C  3(C  Jp  »p  Jp  /p  Jp  /p  *  <p  Jp  »p  <p  3p  >p  <p  /p  <p  O'  ■p  »F  *  ^p 

A=0. 

call  INTEG<CDM,A,B, JNT,Err> 

'  DM=cabs(CDM) 


^  ^  ^  ^  d>  tb  tb  Oj  tb  ^  ^  ^  tb  ^  ^  tb  ^  ^  tb  pb  tb  tb  \b  tb  J/  |b  b  tb  ^  ^  ^ 

^  ^  ^  ^  ^  ^  #P  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  O'  ^  *  *p  0*  *p  ^  ^  ^  ^ 

He  Perfect  dielectric  method  # 

*  E  fields  <  *r  )  * 

*  E<1)  theta  pol.  * 

t  E<2)  horizontal  pol ♦  # 

da  ^  ^  ^  tb  ^  ^  tb  ^  ^  tb  lb  lb  O/  ^  ^  tb  'L'  tb  J/  »X>  tb  ^  tb  <b  ^  lb  tb  tb  tb  ^ 

Ol  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^P  ^  ^  ^  <p  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *P  ^ 

E  <  1 )  =30  a  #Ko*TT*DM*c os  <  THg  >  #cos  ( PHI  r ) 
E<2)=30.*Ko*TH*DM*sin<PHIr> 

Upd=l/<2.#ETAo)*E<  J)*E< J)#C0S<THa)/C0S<THg >*<  (Na/Ng  )**2> 
Ubs=Upd*cos  <  THa ) /cos  <  THg  > 


return 


SUBROUTINE  INTEG< INTGRL, A, B, ARG,Errg ) 

C  INTEGRATION  SUBROUTINE 

^  df  ^  ^  Uj  dr  ^  ^  ^  dr  d'^^^vk\h-^dr\bibdr^^^^dr^d>  dr  \As  d-  dr  ^  ^  dr  dr  ^  ^  dr  tb  »y  ^ 

r^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  rp  ^  t*  ^  ^  ^  n**  rp  ^  rp  ^  *  n*  ^  ^  ^  ^  ^  ^  ^  ^  ^p  *p 


#  Input  # 

#  A,B  endpoints 

#  ARG  integrand  function  # 

#  Errg  greatest  integration  error  * 

#  Output  # 

#  INTGRL  value  of  integral  * 

#  Errg  * 


#*  #  #  *  *#  #**  #  He  *  *  *)K  #  *  *  #  ##  *»oK  #*  *  *:  #  ;f:  ##  #  ft##  *  *  **#*  * 
EXTERNAL  ARG 

COMPLEX  CUM, DEL , Y (5>  , ARG, INTGRL ,YIU 
DOUBLE  PRECISION  RDP,XN,RANG 

RANG=B-A 
NINIT=40 
NSTEPS=N.TNIT 
ERR0R=0.0 
ERRMIJL=1 .0/180.0 
DO  3  1=1,5 
3  Y ( I ) =0.0 

CUM=0.0 
XN=2*NSTEFS 
RDP-A  ' 

R=RDP 

Y  < 1 )=ARG<R) 

15  DO  20  1=2,5 

RDP=RDP+RANG/XN 

R=RDP 

Y  < I )=ARG<R) 

20  CONTINUE 

25  YIV=Y<1 >+Y<5)-4.0*<Y<2)+Y<4>  >+6.0*Y<3> 
ERR0R=ERRMUL*CABS ( YI V ) 

26  IF(ERROR.GT.Errg)  Errg=ERR0R 

DEL=( Y  <l)+Y<5)+4.0*(Y(2)+Y<4) >+2.0#Y (3) > 
DEL=RANG*DEL/ <  3 . *XN ) 

CUM=CUM+DEL 

IF<<R+1.0F-5).GT.B)  GO  TO  80 
Y<1)=Y<5) 

DO  30  1=2,5 

RDP=RDP+RANG/XN 

R=RDP 

Y  < I )=ARG<R) 

30  CONTINUE 

GO  TO  25 

80  INTGRL=CUM 
RETURN 
end 


COMPLEX  FUNCTION  INT<S> 

*  INTEGRAND  FUNCTION 

|i 

^  ^  d#  jb  ^  ^  ^  Uj  ^  ^  dl  U>  tV  ^  ^  iV  lb  ^  ^  ^  ^  d-  ^  *Xr  ^  ^  ^  ^  ib  di  ^  ^  ^  ^  ^  ^  ^  lb  ^  ^  lb  ^  lb  d>  lb  ib  ib  ^  >b  ib  ^  tb 

r^  1^  ^n  *T^  1^*  ^P  ^P  ^n  ^P  ^P  ^P  ^P  '  *P  ^P  ^P  ^P  ^P  ^P  ^P  ^P  ^P  ^P  ^P 


#  Input  from  INTEG  * 

t  S  distance  on  antenna  * 

#  Input  from  MAIN  * 

#  Kr  relative  propagation  constant  of  antenna  * 

#  ALP  attenuation  constant  of  antenna  * 

#  shift  phase  shift  for  second  half  of  array  # 

#  Input  from  method  subroutine^)  # 

#  THa  angle  from  the' normal  in  air  < radians)  # 

#  PHIr  azimuthal  angle  (radians)  * 


ib  ^  ^  ^  ^  ^  d/  ib  \ V  vb  ib  d/  ^  ib  df  ib  ib  «b  d/  d/  ^  d/  ibtbtb  ib  >b  d>  ib  ^  ib  tb 

^  ^  ^  ^  ^  ®  ^  »r  •p  *p  <p  *p  n*  ^  *p  *•  *•  rp  A*  ^  ^  ^  #p  An  ^  *  *  An  ^  ^  »p  *  ^  ^  ^  ^  Ah  ^  ^  ^  ^  ^  »p  ^  ^  ^  ^  ^  <p 

COMMON  /b 1/ ALP, Kr, shift 
COMMON  /b2/THa ,PHI r 
REAL  Kr ,LAMo 
COMPLEX  ZJ,I 


y  ^  ^  ^  j|^  y  |j^  ^  lj^  <|^  lj^  lj^  ^  j|^  ^  ^  <|^  ^  ^  lj^  ^  \jf 

*  FQ  operating  frequency  * 

*  C  speed  of  light  # 

^  ^  ^  ^  ^  ^  ^  ^  ,p  *  ^P  ^  A'  ^  ^  ^  ^  ^  rp  ^  )p  ^  ^  A’  *P  ^  ^  ^ 

data  FG,C,ZJ, PI/37. 5e6,3e8, CO. , 1 ♦ ) ,3. 1415926536/ 

dr  ^  <b  ^  ^  ^  ^  ^  lb  ^  tb  ^  ib  ^  ^  d>  \b  ^  ^  ^  ^  ib  ^dfib  d,  lb  U/  \L*  lb  ^  ^  lb  ^  dr  ^  dr  ^  ^  ^  ib  ^  ib 

*  LAMo  -  free  space  wavelength  # 

*  BET  -  antenna  propagation  constant  * 

*  BETo  -  free  space  propagation  constant  * 

d(  d/  dr  ^  d#  ^  ^  ^  ^  lb  ^  lb  «W  tb  ^  vV  ib  tb  ^  tb  ^  ib  lb  ib^ib  lb>  ib  ^  U>  ib  -b  lb  lb  dr  ^  ib  ib  tb  ^  ib  ^ 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *p  ^  ^  rp  *p  ^  ^  rp  ^  ^  ^  ^  ^  /p  ^  ^  ^  ^  ^  ^  ^  An  rp  /p  i^\  ^  ^  ^  ^  ^  .y.  Ah 

LAMo=C/FQ 

BET=Kr*2.#PJ/LAMo 

BETo=2.*PI/LAMo 

lb  d,  ^  d*  ^  ^  ^  d,  ^  ib  ^  ^  lb  ^  ^  lb  ^  ^  ^  ^  ^  dr  lb  ^  ib  ^  ^  ib  ib  ^  ^  ib  ^  ^  ib  )b  ^ 

*  ^  ^  ^  ^  ^  ^  rp  *  ^  rp  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *  »p  ^  ^  rp  ip  ^  ^  ^  ^  ^  ^  ip  ^  ^ 

*  special  instructions  for  the  array  * 

^  ^  An  ^  ^  ^  ^  ^  ni  ^  A"  An  ^p  ^  ^n  ^n  ^  *  ^p  *  ^  ^  ^  ^  ^  ^p  ^p  ^p  ^p  ip  ^  ^P  An  ^n  ^P  ^  * 

ifCS.gt.16.)  then 

ifCS.lt. 20.)  goto  10 
d=36-S 

gam=shift*PI/180. 
else 
d=S 
gam=0. 
end  if 


t  *  ft  *  *  He  ft  ft  ft  *  *  ft  *  *  *  *  *  *  ft  ft  ft  ft  ft  ft  ft  *  * 

*  I  -  CURRENT  ON  ANTENNA  * 

*ftftft**ftft#*ftftft*ftftftft#ftft#*ftftft# 

I  =SI  N  C  BETftd  )  ftCEXP  <  -  Ai.  P*(16.-d)  +7 Jftgam ) 


-»  ->  -»  -■ 


B-6 


*  INT  -  DIPOLE  MOMENT  # 

*****************  ********  ****** 

INT=I*CEXP  <  7 J*BETo#S*SIN  < THa ) *COS  <  PHI r ) ) 

RETURN 


10  INT=0 
return 
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